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Abstract. We have collected one-dimensional raster-scan observations of the active star-forming region Sharpless 
171 (S171), a typical Hn region-molecular cloud complex, with the three spectrometers (LWS, SWS, and PHT-S) 
on board ISO. We have detected 8 far-infrared fine-structure lines, [Om] 52 /im, [Niii] 57 /im, [O i] 63 (jm, [Om] 
88 /im, [N n] 122 fim, [O i] 146 fj,m, [C n] 158 /im, and [Si n] 35 /im together with the far-infrared continuum and the 
H2 pure rotation transition (J = 5-3) line at 9.66 /im. The physical properties of each of the three phases detected, 
highly-ionized, lowly-ionized and neutral, are investigated through the far-infrared line and continuum emission. 
Toward the molecular region, strong [0 1] 146 fim emission was observed and the [O 1] 63 fim to 146 fim line ratio was 
found to be too small (~ 5) compared to the values predicted by current photodissociation region (PDR) models. 
We examine possible mechanisms to account for the small line ratio and conclude that the absorption of the [0 1] 
63 /jm and the [C 11] 158 (im emission by overlapping PDRs along the line of sight can account for the observations 
and that the [0 1] 146 /im emission is the best diagnostic line for PDRs. We propose a method to estimate the 
effect of overlapping clouds using the far-infrared continuum intensity and derive the physical properties of the 
PDR. The [Si 11] 35 /im emission is quite strong at almost all the observed positions. The correlation with [N 11] 
122 (im suggests that the [Sill] emission originates mostly from the ionized gas. The [Sill] 35 /im to [Nil] 122 /im 
ratio indicates that silicon of 30% of the solar abundance must be in the diffuse ionized gas, suggesting that 
efficient dust destruction is undergoing in the ionized region. 
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1. Introduction 

Far- infrared (FIR) spectroscopy of the interstellar medium 
(ISM) provides us a great deal of information on the na- 
ture of the ISM. For Hn regions, FIR forbidden lines are 
useful tools to investigate the physical properties, such as 
the electron density and elemental abundance. They are 
less subject to extinction and less sensitive to the elec- 
tron temperature than optical forbidden lines. Rubin et 
al. 1)1994[1 described a semiempirical methodology to de- 
rive the electron density, the effective temperature for the 
ionizing star, and the gas-phase heavy element abundance 
from FIR line emissions based on the ionization bounded 
models. In neutral regions, in which the gas is warmer than 
the molecular gas observed in radio frequencies, the kinetic 
temperature of atoms is high enough to excite FIR emis- 
sion and a number of forbidden lines are emitted. In the 
interface region between the ionized and molecular gas, in- 
tense far- ultraviolet (FUV)(6 eV < hv < 13.6 eV) photons 
photodissociate molecules and photoionize heavy elements 
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with ionization potential less than the Lyman limit. This 
region is called the photodissociation region (PDR) , where 
most energy is emitted in the FIR lines and continuum. 
Theoretical models have been investigated for PDRs with 
various physical conditions (Tielens & Hollenbach 1985 
Hollenbach et al. lTMD Wolfire et al. HMfl Kaufman et al. 
ITMfl see Hollenbach & Tielens for a review). These 
models take account of the chemistry and the energy bal- 
ance and predict the cooling line intensities from PDRs. 

In this paper, we report the results of a spectroscopic 
investigation of Sharpless 171 (S171) with the Infrared 
Space Observatory (ISO; Kessler et al. I1996[) . S171 is a 
large Hn region and is the main cloud associated with 
the Cepheus OB4 stellar association (Yang & Fukui l992). 
The star cluster Be 59 located at the central portion of the 
nebulous region is the ionizing source of this region. It con- 
tains one 07 star and several later-type stars. The phys- 
ical properties of the ionized gas have been investigated 
in radio continuum and recombination lines (e.g., Felli et 
al. 115771 Rossano et al. 115501 Harten et al. HMl} . Yang 
& Fukui l|1992[) have mapped the large-scale molecular 
gas distribution with the 13 CO (J = 1-0) line, suggesting 
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Fig. 1. 13 CO (J = 1-0) integrated intensity map (in gray- 
scale) and the 0.61 GHz radio continuum emission map (in 
contours) (Yang & Fukui 1992). The position of Be 59 is 
shown by the plus (+). The thick solid line indicates the 
line along which the present observations were carried out. 

that Be 59 generates the ionization front on the surface of 
two dense molecular clumps (CI and C2). They suggested 
from the dynamics of the molecular clumps that the dense 
gas is contacting with the continuum source and the ion- 
ization front is driving shocks into the CI clump. S171 is 
a typical transition region from ionized gas to molecular 
clouds. It has a scale of several tens arcminutes, which is 
appropriate for mapping observations in the FIR to study 
the physical properties of the ionized gas and PDR com- 
plex. 

In Sect. 2, the observations and the data reduction 
are described. The results are presented in Sect. 3 and 
discussed in Sect. 4. A summary is given in Sect. 5. 

2. Observations and data reduction 

We made one-dimensional raster-scan observations of 
the S171 region with the Long- Wavelength Spectrometer 
(LWS; Clegg et al. HMty . the Short- Wavelength 
Spectrometer (SWS; de Graauw et al. I1996f> . and the 
PHT Spectrometer (PHT-S; Lemke et al. Laureijs 
et al. I2UU2|) on board ISO. The one-dimensional raster- 
scan started from the position near Be 59 to the region 
of ~ 5 pc from Be 59, which is spanned from the ion- 
ized gas to the molecular cloud (Fig. ^) . The observations 
were carried out at 24 positions (Table GJ, each of which 
was separated by 50". The full grating scan mode AOT 
LWS01 was used in the LWS observations to cover the 
wavelength range from 43 to 197 /jm with A/AA ~ 100- 
300 (Gry et al. E)02). Four grating scans were carried 
out for each raster position with the spectral sampling 
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Fig. 2. LWS spectra of the 24 raster positions. The spectra 
are shifted in the y-axis for clarity from top (position 1) 
to bottom (position 24). One division mark in the y-axis 
corresponds to 10~ 6 Wm -2 ^m _1 sr~ 1 . 

of a half of the resolution clement. The total integration 
time at a grating position was 2 sec. The LWS beam size 
was 66"-86", depending on the wavelength, which corre- 
sponded to 0.27 ~ 0.35 pc at the distance of S171 of 850 pc 
(MacConnell 1968), and the beam overlapped with each 
other. The Off-Line Processing (OLP) version 10.1 data 
obtained from the ISO Archival Data Center were used for 
the present study. The spectra were further processed by 
using the ISO Spectral Analysis Package (ISAP 1 ). After 
removing remaining noises, the spectrum at each posi- 
tion was averaged over individual spectral scans. The long 
wavelength part of the LWS spectra of extended sources 
is affected by the presence of fringes in the continuum 
(Swinyard et al. I1998|l . therefore for A > 90 fim the spec- 
tra have been corrected for fringes using the procedure 
within ISAP. Finally, the extended source correction was 
applied and the fluxes were converted into surface bright- 
ness according to the beam size of the detectors. We used 

The ISO Spectral Analysis Package (ISAP) is a joint de- 
velopment by the LWS and SWS Instrument Teams and Data 
Centers. Contributing institutes are CESR, IAS, IPAC, MPE, 
RAL and SRON. 



Table 1. Summary of the line intensity of the LWS observations 



position line intensities (10 s Wm 1 sr r j°' 

No. I (degree) 6 (degree) d (pc) fc [Om]52jtm [Nm] 57^m c [Ol]63^m [O m] 88 /mi [Nil] 122 /im [O i] 146 /tm d [Cll]158/mr 



1 


118 


.225 


4 


.987 





.14 


31.50 


± 


2.88 


10.73 


± 


1.26 


< 1.01 


52.82 


± 





.65 





.73 


± 





.17 




< 


0.24 


6.49 ± 0.26 


2 


118. 


.211 


4 


.987 





.15 


35.51 


± 


3.67 


13.17 


± 


1.23 


< 1.16 


57.01 


± 


1 


.63 





.84 


± 


0. 


.23 




< 


0.20 


7.37 ± 0.26 


3 


118 


.197 


4 


.986 





.33 


32.90 


± 


4.84 


14.25 


± 


1.13 


2.84 ± 0.55 


59.16 


± 


1 


.13 


1 


.23 


± 


0. 


.19 





.31 


± 0.16 


9.44 ± 0.36 


4 


118. 


.183 


4 


.986 





.53 


32.52 


± 


4.72 


12.85 


± 


1.18 


8.49 ± 0.56 


58.48 


± 


1 


.11 


1 


.31 


± 


0. 


.24 





.50 


± 0.13 


15.82 


± 0.61 


5 


118. 


.169 


4 


.986 





.74 


43.68 


± 


2.93 


15.99 


± 


1.17 


15.93 ± 0.63 


65.29 


± 


1 


.15 


1 


.39 


± 


0. 


.36 





.97 


± 0.18 


20.86 


± 0.88 


6 


118. 


.155 


4 


.986 





.94 


47.48 


± 


3.89 


18.17 


± 


1.18 


20.66 ± 1.20 


75.19 


± 


1 


.66 


1 


.40 


± 





.19 


1 


.36 


± 0.16 


23.99 


± 0.65 


7 


118 


.141 


4 


.986 


1 


.15 


43.30 


± 


5.03 


17.64 


± 


0.79 


13.64 ± 1.25 


75.05 


± 


1 


.96 


1 


.73 


± 


0. 


.32 


1 


.10 


± 0.13 


24.13 


± 0.82 


8 


118. 


.127 


4 


.985 


1. 


.35 


38.12 


± 


3.88 


16.35 


± 


0.96 


9.92 ± 0.86 


67.45 


± 


1 


.65 


1 


.99 


± 


0. 


.37 





.80 


± 0.12 


24.20 


± 0.97 


9 


118. 


.114 


4 


.985 


1. 


.56 


38.10 


± 


3.41 


14.97 


± 


1.02 


7.83 ± 0.61 


61.18 


± 


1 


.72 


2 


.20 


± 


0. 


.24 





.65 


± 0.12 


23.72 


± 0.87 


10 


118. 


.100 


4 


.985 


1 


.76 


34.52 


± 


3.88 


13.82 


± 


1.04 


7.82 ± 0.38 


59.02 


± 


1 


.87 


1 


.84 


± 


0. 


.23 





.49 


± 0.13 


23.80 


± 0.99 


11 


118. 


.086 


4. 


.985 


1. 


.97 


31.66 


± 


3.44 


12.05 


± 


1.00 


7.16 ± 0.58 


49.63 


± 


1 


.92 


1 


.54 


± 


0. 


.23 





.56 


± 0.10 


24.28 


± 1.08 


12 


118. 


.072 


4 


.984 


2. 


.17 


30.08 


± 


1.52 


12.36 


± 


0.69 


6.54 ± 0.41 


46.47 


± 


2 


.75 


1 


.55 


± 





.24 




< 


0.30 


21.88 


± 0.95 


13 


118. 


.058 


4 


.984 


2 


.38 


29.73 


± 


2.35 


10.48 


± 


0.66 


5.32 ± 0.46 


45.13 


± 


2 


.57 


1 


.50 


± 





.23 





.35 


± 0.10 


20.86 


± 1.02 


14 


118. 


.044 


4 


.984 


2 


.58 


22.66 


± 


2.08 


10.66 


± 


0.62 


6.61 ± 0.55 


39.52 


± 


2 


.15 


1 


.75 


± 





.27 





.39 


± 0.09 


23.01 


± 1.05 


15 


118 


.030 


4 


.984 


2 


.79 


17.73 


± 


1.46 


9.09 ± 0.66 


8.20 ± 0.47 


34.19 


± 


1 


.89 


1 


.46 


± 


0. 


.27 





.52 


± 0.14 


25.23 


± 1.25 


16 


118 


.016 


4 


.984 


3 


.00 


26.40 


± 


3.35 


8.71 ± 0.66 


10.23 ± 0.44 


34.02 


± 


1 


.82 


1 


.17 


± 


0. 


.29 





.74 


± 0.13 


26.94 


± 0.96 


17 


118. 


.002 


4 


.983 


3 


.20 


18.74 


± 


1.45 


9.29 ± 0.56 


13.13 ± 0.52 


31.12 


± 


1 


.79 


2 


.18 


± 


0. 


.31 





.95 


± 0.12 


32.12 


± 1.40 


18 


117 


.988 


4 


.983 


3 


.41 


17.73 


± 


1.49 


8.47 ± 0.55 


17.23 ± 0.51 


27.12 


± 


1 


.82 


2 


.92 


± 





.45 


1 


.64 


± 0.17 


34.89 


± 1.57 


19 


117. 


.974 


4 


.983 


3. 


.61 


14.49 


± 


2.20 


6.84 ± 0.86 


19.56 ± 1.33 


23.29 


± 


1 


.33 


2 


.44 


± 





.49 


3 


.28 


± 0.20 


40.52 


± 1.60 


20 


117 


.960 


4 


.983 


3 


.82 


10.50 


± 


1.97 


5.60 


± 0.77 


16.30 ± 0.77 


19.70 


± 


1 


.45 


2 


.15 


± 


0. 


.39 


3 


.57 


± 0.29 


37.33 


± 1.51 


21 


117. 


.946 


4 


.982 


4. 


.03 


9.67 ± 2.28 


5.01 ± 0.60 


17.68 ± 0.66 


15.36 


± 


1 


.33 


1 


.91 


± 


0. 


.54 


4 


.92 


± 0.26 


37.97 


± 1.20 


22 


117 


.932 


4 


.982 


4. 


.23 


6.86 ± 1.70 


< 


1.70 


14.46 ± 0.93 


12.01 


± 





.86 


1 


.01 


± 





.32 


3 


.44 


± 0.33 


32.48 


± 1.97 


23 


117. 


.918 


4 


.982 


4, 


.44 


6.19 ± 1.46 


2.50 ± 0.51 


8.58 ± 0.67 


10.09 


± 





.85 





.73 


± 





.16 





.82 


± 0.11 


22.00 


± 1.11 


24 


117. 


.904 


4 


.982 


4. 


.64 


3.15 ± 1.30 


3.25 ± 0.57 


6.35 ± 0.57 


9.23 ± 1.04 





.85 


± 


0. 


16 





.79 


± 0.25 


17.21 


± 0.75 



a Upper limit (3a) is estimated from the statistical error of the baseline. 
b Distance from Be 59. 

c Weighted mean value of the intensities determined from the SW2 and SW3 channels. 
d Weighted mean value of the intensities determined from the LW3 and LW4 channels. 
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Table 2. Summary of the line intensity of the SWS and PHT-S observations 



No. 


[Si II] 35 Aim 


line intensities (10 
H 2 9.66 /mi (SWS) 


-"Wm-'sr 1 )" 

H 2 9.66 Mm (PHT-S) 


H2 6.9 /Urn 


1 


< 4.37 


< 3.02 


< 


3.94 


< 2.52 


2 


< 6.08 


< 2.20 


< 


4.88 


< 5.62 


3 


< 6.83 


< 2.22 


< 


4.53 


< 3.80 


4 


< 7.10 


< 2.44 


< 


5.30 


< 3.82 


5 


< 4.18 


3.05 ± 1.03 


< 


6.16 


< 7.07 


6 


7.03 ± 1.96 


< 2.74 


< 


6.72 


< 6.08 


7 


6.20 ± 1.68 


< 2.28 


< 


6.80 


< 6.08 


8 


5.55 ± 1.96 


< 2.19 


< 


4.37 


< 5.81 


9 


7.24 ± 1.40 


< 2.29 


< 


3.51 


< 4.58 


10 


5.73 ± 1.81 


< 2.23 


< 


3.61 


< 5.63 


11 


7.25 ± 1.99 


< 2.55 


< 


4.71 


< 5.62 


12 


< 6.07 


< 2.26 


< 


4.88 


< 4.89 


13 


< 6.75 


< 2.27 


< 


3.85 


< 5.28 


14 


< 5.33 


< 2.00 


< 


5.90 


< 7.04 


15 


6.47 ± 1.82 


< 3.64 


< 


6.01 


< 5.98 


16 


< 5.76 


< 1.81 


< 


4.45 


< 5.62 


17 


7.97 ± 1.72 


< 2.64 


< 


5.30 


< 6.34 


18 


10.62 ± 1.71 


< 3.10 


< 


4.19 


< 8.70 


19 


17.76 ± 1.92 


2.81 ± 0.87 


< 


6.76 


< 14.13 


20 


11.05 ± 1.65 


5.82 ± 0.63 


< 


6.93 


< 17.39 


21 


18.62 ± 2.20 


8.34 ± 1.09 


< 


6.45 


< 21.74 


22 


7.86 ± 1.79 


3.30 ± 0.97 


< 


6.84 


< 14.31 


23 


< 4.18 


< 3.87 


< 


4.45 


< 7.97 


24 


3.87 ± 1.41 


< 2.37 


< 


5.56 


< 7.61 



a Upper limit (3<r) is estimated from the statistical error of the baseline for the SWS and from the uncertainty in the flux of 
the corresponding detector channel for the PHT-S 



the latest values of the extended source correction factors 
and the beam sizes given in Gry et al. I|2002[l . The line 
intensities were obtained by Gaussian fits. The absolute 
flux calibration uncertainty is reported to be 10-20% for 
point sources and 50% for extended sources (Gry et al. 
2002]. The relative accuracy is assumed to be 30%, which 
was estimated from the gaps in the signal levels of the 
detectors. 

The same positions were observed with the SWS and 
the PHT-S. For the SWS observations the grating line 
profile mode AOT SWS02 was used to observe [Sill] 
35 /im and the H2 pure rotational transition (J = 5-3) 
at 9.66 /mi. The aperture size was 20" x 33" and 14" x 20" 
and the spectral resolution was about 1500 and 2000 for 
the [Sin] 35 /im and H2 9.66 /im lines, respectively (Leech 
et al. I2002|l . The OLP 10.1 data were obtained from the 
ISO Archival Data Center and used in the present study 
except for the data at the position 23 in Tabled for which 
only the processed data of version 9.1 were available. At 
each raster position one up-down grating scan was carried 
out with the total integration time of 200 sec. We used the 
ISAP for further data processing and obtained the line in- 
tensities by Gaussian fits. The flux calibration accuracy is 
reported to be 7% and 22% for the H 2 9.66 fim and [Sill] 
35 /im lines, respectively, for point-like sources. The con- 
version factors for extended sources are estimated to be 
accurate within 10% (Leech et al. [2002). 



The PHT-S had a resolving power of A/AA ~ 90 in 
two wavelength bands with the aperture size of 24" x 24". 
The data were reduced by the Phot-Interactive Analysis 
(PIA 2 ) version 8.0. The total integration time was 128 sec 
at each raster position. The absolute flux calibration accu- 
racy is about 10% or less for point sources (Laureijs et al. 
2002). Because of the low spectral resolution, we could not 
apply Gaussian fits to estimate the line intensity. Instead, 
we simply measured the flux of the detector channel at 
the line center and subtracted the continuum estimated 
from the adjacent detectors. 

3. Results 

Figures^ andHshow the LWS, SWS, and PHT-S spec- 
tra of S171, respectively. We detected 7 forbidden lines, 
[Om] 52 /im, [Nm] 57 /im, [Oi] 63 /im, [Om] 88 /im, [Nil] 
122 /im, [Oi] 146 /im, and [C11] 158 /im in the LWS spec- 
tra of almost all the observed positions. The results of the 
LWS observations are summarized in Tabled The [Nm] 
57 /im and [Oi] 146 /^m lines were observed by two adja- 

2 The ISOPHOT data presented in this paper were re- 
duced using PIA, which is a joint development by the ESA 
Astrophysics Division and the ISOPHOT Consortium with the 
collaboration of the Infrared Processing and Analysis Center 
(IPAC). Contributing ISOPHOT Consortium institutes are 
DIAS, RAL, AIP, MPIK, and MPIA. 
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Fig. 3. SWS spectra of the 24 raster positions. The spectra 
are shifted in the //-axis for clarity from top (position 1) 
to bottom (position 24). One division mark corresponds 
to 10~ 5 Wm~ 2 /mi^sr -1 . 



cent detectors, SW2 and SW3, and LW3 and LW4, respec- 
tively. The line intensities derived from different detectors 
are in agreement with each other within the estimated 
uncertainties. We use the weighted mean of the intensities 
from the two detectors in the following analysis. In the 
case that the line was detected only by one detector with 
the other being an upper limit, we adopted the intensity 
of the detected channel. [O in] 88 /im was also observed by 
two adjacent detectors, SW5 and LW1. We used only the 
SW5 data for [O in] 88 /im because of its higher spectral 
resolution. We have detected [Sin] 35 /im and the H2 pure 
rotational transition at 9.66 /im in the SWS observations, 
while H2 9.66 /im and 6.9 /im (J = 7-5) were not detected 
at the 3a level in any PHT-S spectra. Upper limits for 
these lines in the PHT-S spectra were estimated from the 
uncertainties in the flux of the corresponding detectors. 
The upper limits of the 9.66 /im line are compatible with 
the SWS results except for position 21, where the discrep- 
ancy may be attributed to the difference in the aperture 
size if the emission is highly peaked. The results of the 




Wavelength [urn] 



Fig. 4. PHT-S spectra (5.6-11.6/im) of the 24 raster po- 
sitions. The spectra are shifted in the j/-axis for clarity 
from top (position 1) to bottom (position 24). One divi- 
sion mark corresponds to 10~ 6 Wm~ 2 /im _1 sr _1 . 



SWS and PHT-S observations are given in Table [3 The 
errors of the line intensities include both the fitting er- 
rors and the statistical errors of the baseline for the LWS 
and SWS data. Uncertainties in the absolute flux are not 
included in the errors in Tables ^ and [21 For the line emis- 
sion less than 3 times the uncertainty (3<t) we regard it as 
non-detection and give 3er as the upper limit. 

Figure shows the spatial distributions of the line in- 
tensities as a function of the distance from Be 59, d. The 
[Oi] lines and the [Nil] 122 /im line show spatial distribu- 
tions similar to the 0.61GHz radio continuum flux (Fig.[6j). 
The intensities of these lines show maxima at d ~ 1 pc and 
~ 4 pc, indicating the presence of high density gas at these 
positions. The [Cli] 158 /im, [Sin] 35 /im, and H2 9.66 /im 
emissions also show a strong peak at ~ 4 pc, but their 
intensity distributions become rather flat around ~ 1 pc. 
The [O ill] and [N ill] lines appear strong at d ~ 1 pc and 
monotonically decrease with d, indicating a decrease in 
the abundance of the doubly ionized ions. 

The coexistence of ionic lines of different ionization 
degrees as well as neutral atomic and molecular lines in- 
dicates that the observed region contains gases of various 
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Fig. 5. Observed line intensities against the distance from Be 59. The arrows indicate the 3<r upper limit, a) [Oin] 
52 /xm (lower points of +) and 88 (im (upper points of *); b) [Nm] 57 /im; c) [Oi] 63 /xm (upper points of +) and [Oi] 
146 /xm multiplied by 2 (lower points of *); d) [Nil] 122 /xm; e) [Cn] 158 /xm; f) [Sill] 35 /xm; g) H 2 9.66 /xm by the 
SWS; h) H 2 9.66 ^m by the PHT-S; and i) H 2 6.9 /xm. 



physical conditions. In the following we attribute each ob- 
served line to three distinct physical phases, the highly- 
ionized gas, the lowly-ionized gas and the PDR gas, and 
investigate the properties of each gas. 

3.1. Highly-ionized gas 

[Oni] 52 /xm, 88 /xm and [Nm] 57 /xm originate from the 
highly-ionized gas. The [Oni] 52 /xm to 88 /xm intensity 
ratio is in general a function of the electron density and 
temperature (Osterbrock 1989). However, this ratio is 
insensitive to the electron temperature around 10 4 K, 
and the electron density can be estimated directly from 
the observed ratio. The maximum extinction in the ob- 
served positions is estimated as T cxt ~ 0.1 at 52 /xm 
based on the 12 CO line intensity at the CI cloud (Yang 
& Fukui I1992fl and the standard dust model (Draine & 
Lee [TTh 1 i with the conversion factor of 2.8 x 10 20 IT 
molecule K _1 km - s cm" " 2 (Mizutani et al. ESUUl for de- 
tails). Thus the extinction by dust does not affect the ob- 
served 52 /xm to 88 /xm ratio. 
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Figure [7| shows that the observed ratio is approxi- 
mately ~ 0.6 except for position 24 where the detection 
of [O in] 52 /jm is marginal. The mean ratio 0.6 indicates 
n e ~ 20 cm -3 for T e = 10 4 K with the parameters for 
the transition given in Osterbrock 1|1989[1 . The low den- 
sity limit of the ratio is 0.56 for T e = 10 4 K. The uncer- 
tainty in the collision coefficients leads to an uncertainty of 
about 20% in the ratio at the low density limit (Mizutani 
et al. l2002[) . The relative flux calibration error leads to an 
uncertainty of ±0.2 in the ratio. Taking account of these 
uncertainties, the upper limit at the electron density is es- 
timated as 150 cm" 3 . The lower limit is not constrained. 
The derived range of the electron density (n e < 150 cm~ 3 ) 
is about the same as that in the Carina region (n e < 100 
cm~ 3 ; Mizutani et al. 2002), indicating the presence of the 
diffuse Hn region also in S171. According to Mizutani et 
al. H2002f) . the line-of-sight thickness of the line-emitting 
gas L can be given by 



1.2 x 



J([0 in]88/zm) 



20 cnT 



pc, (1) 



0.5 



when n e is well below the critical density, the fill- 
ing factor is unity, and log(0/H) = —3.3 (Holweger 
2001). For the maximum [O in] 88 /im intensity of ~ 
7.5 x 10~ 7 Wm~ 2 sr _1 , L becomes approximately 9 x 
(20cm~ 3 /n e ) 2 pc. This thickness is comparable with the 
projected size of S171 of several pc, suggesting that the 
highly-ionized gas is associated with S171. 

The [N in] 57 /im to [O ill] 52 /im intensity ratio is also 
insensitive to the temperature, and thus is a function of 
the electron density and the abundance ratio of N ++ and 
++ (Mizutani et al. I2002[) . N and O have similar sec- 
ond ionization potential energies and both lines can be 
assumed to come from the same region. From the ob- 
served ratio of ~ 0.4 at most positions, we derive the 
abundance ratio N++/0++ as 0.23 or 0.29 for T e = 10 4 K 
and n e — 20 cm~ 3 or n e = 150 cm~ 3 , respectively. Taking 
account of the relative flux calibration error of the LWS, 
we estimate N++/0++ = 0.2 ± 0.1 for n e < 150 cm~ 3 . 
This is in agreement with the value of about 0.2 at the 
galactocentric distance i? ga i — 8.5 kpc, which is estimated 
from 34 compact Hn regions located between i? ga i = 0.3 
and 15 kpc (Martin-Hernandez et al. 20021), and is also 
similar to the value of ~ 0.19 of the core ionized gas in 
the Carina nebula (Mizutani et al. [2002 ) . 

EE 



According to Rubin ct al. (1994) the effective temper- 
ature of the ionizing star can be estimated from the [N III] 
57 /im to [N II] 122 /im line ratio for the ionization bounded 
condition. With T e = 10 4 K and n e < 150 cm -3 , we de- 
rive (N ++ )/(N + ) ~ 0.69-0.86 from the total intensities of 
[Nmj 57 /im and [Nil] 122 /mi of the 24 positions. This 
indicates T e g ~ 37000-37500 K of the ionization source, 
which is cooler of about 1000 K than that of the earliest 
type star in the S171 region, 07 (T off = 38500 K). The 
[Nil] 122 /im emission comes also from the diffuse ion- 
ized gas and the ionization bounded condition may not 
hold. Conversely an 07-type star can provide a sufficient 
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52 /xm to 88 /im ratio against the distance 



amount of the [N in] 57 /xm emission to account for the ob- 
served intensity indicated by the [Nil] 122 /im emission. 

3.2. Neutral region 

The [Oi] 63 /im, 146 /im, [Cn] 158 /im, [Sin] 35 /im, and 
H2 lines originate from the PDR, where the FUV flux is 
converted into atomic and molecular gaseous line emis- 
sions by the photoelectric heating. We estimate the physi- 
cal properties from the observed intensities of the line and 
continuum emission, using the PDR model by Kaufman et 
al. (1999). The major model parameters are the gas den- 
sity, n, and the radiation field strength, Go, in units of the 
solar neighborhood (1.6 x 10" 6 Wm" 2 ; Habing tlSSgj l. The 
model gives integrated line intensities seen in the face-on 
view. 

First, we fit the continuum emission of the LWS spec- 
tra with a gray-body radiation as 



F(X) = tioo x 



A 



100/im 



-0 



Bx(T d ), 



(2) 



where Td is the dust temperature, tioo is the optical depth 
at 100 /an, B\ is the Planck function, and j3 is the power- 
law index of the emissivity on the wavelength. We assume 
(3 = 1 because it gives the best fit. The obtained Td ranges 
from 34 K to 53 K and tends to decrease with the distance 
from Be 59. The total far-infrared flux FIR is given by 



FIR 



F(X) dX 



A 



100/im 



Bx(T d ) dX. 



(3) 
(4) 



In practice, the integration was carried out over the spec- 
tral range that makes effective contributions. The far- 
infrared flux should be balanced by the absorbed energy 
in the ultraviolet to optical region as 



G uv x (1.6 x lO^Wm^) 



47rj;dA 



(5) 



912.1 
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Fig. 8. Radiation field intensity Guv estimated from the 
luminosity of an 07 type star at the origin (solid line) and 
that from the temperature by fitting the LWS spectrum 
(crosses with errorbars, see text). 

= -j^FIR oc 7^+ 4 , (6) 

(Tabs) 

where is the mean intensity from the heating sources 
and (r a b s ) is the weighted mean of the absorption optical 
depth over J£. We assume (T a b s )/Tioo = 700 (Okumura 
IT333I Chan et al. fMty . In Fig. |H| the crosses with er- 
rorbars show the derived Guv- We also calculate Guv as 
L/iircP, where L (= 2.6 x 10 5 Lq) is the luminosity of 
an 07 type star (solid line in Fig. |HJ. Both estimates are 
roughly in agreement with each other for d > 2 pc. In the 
region closer to the heating sources the estimate from the 
stellar luminosity becomes much larger. It suggests that 
dust grains emitting the FIR emission are not present in 
the vicinity of Be 59. The agreement in the outer region 
supports that Eq. JHJ gives a reasonable estimate for Guv- 
In the PDR model the parameter Go does not include pho- 
tons with the energy of < 6 eV, which do not contribute 
to the photoelectric heating. Thus Go is given by 

,.2066.4 

G x (1.6 x 10" 6 Wnr 2 ) = / „ 4tt J* x d\ . (7) 

J912A 

For an 07 type star with T e g — 38000 K, we estimate 
Guv /Go = 1-3 and use this scaling in the following anal- 
ysis. 

[Cn] 158 /mi and [Oi] 63 fim and 146 /mi are the ma- 
jor cooling lines in the PDR. However [C n] 158 fim comes 
also from the ionized gas. We estimate the density n from 
Go and either of the ratio [O i] 63 fim/FIR (case 1) or the 
ratio [O i] 146 fim/FIR (case 2). With the estimated n, we 
derive the ratio [Cn] 158 fim/ FIR, the surface tempera- 
ture, and the ratio [O i] 146 fim/FIR or [O i] 63 fim/FIR, 
for cases 1 or 2, respectively. The results are plotted in 
Figs. and In both cases the model predicts fairly 
well for d < 3 pc the intensity of the other [O i] line that 
is not used for the input parameter, and always provides 
a smaller intensity of [C n] 158 fim than observed. The ex- 
cess emission of [Cn] 158 fim is considered to come from 



the ionized gas. At d ~ 4 pc, however, the observed [Oi] 
146 fim is much stronger than the model prediction in case 
1 (Fig. Et), or the model predicts too strong [Oi] 63 fim 
in case 2 (Fig. HOh l . The model also predicts the [C n] 
158 fim intensity stronger than observed for d > 3.5 pc 
(case 2; Fig. HOt O or indicates that almost all the [C II] 
158 fim emission comes from the PDR gas for d > 3.5 pc 
(case 1; Fig.|3|i). The [Nil] 122 fim intensity shows a peak 
at d ~ 4 pc and thus there should be a significant con- 
tribution to [Cn] 158 fim from the ionized gas that emits 
[Nil] 122 fim. Therefore the line intensities of [Oi] 63 fim, 
146 fim and [C n] 158 fim together with the FIR continuum 
emission can not be accounted for consistently by the lat- 
est PDR model in either case. In case 2 the gas density has 
a peak at ~ 4 pc (Fig.llOh,). where CO observations indi- 
cate the presence of the main molecular clump CI (Yang 
& Fukui 1992). On the other hand case 1 indicates a local 
minimum in the density and a peak in the temperature at 
~ 4 pc (Figs. & and b), neither of which is compatible 
with the CO observations. Case 2 leads to more consistent 
results with the CO observations (Figs. HHk and b). In the 
following we assume case 2 gives correct results and dis- 
cuss the discrepancy in the [Oi] 63 fim and [Cn] 158 fim 
emission. 

According to the PDR model, [O i] 63 fim is optically 
thick, [Cn] 158 fim is marginally thick (t ~ 1), and [Oi] 
146 fim and the FIR emission are optically thin in a PDR 
(Kaufman et al. 1999). The largest value of the derived 
tioo is ~ 0.005, supporting that the FIR emission is opti- 
cally thin in the PDR of S171. The optical depth effects of 
emission lines are taken into account in the PDR model. If 
several PDR clouds overlap with each other along the line 
of sight, however, the optically thick lines will be attenu- 
ated by foreground PDRs, while optically thin emissions 
will be a simple summation of those from each PDR cloud. 

The FIR intensity is estimated simply as 1.6 x 
10- 6 Guv/4tt (Wm^sr 1 ) in the PDR model by assum- 
ing that all the incident radiation is absorbed and con- 
verted into the FIR emission (Kaufman et al.[1999|. Here 
we define the overlapping factor Z as 



1.6 x 10- 6 (Wm^sr- 1 ) G UV ' 

For Z < 1, Z is supposed to indicate the beam filling fac- 
tor, whereas it should estimate the degree of overlapping 
of clouds for Z > 1 since the FIR emission is optically 
thin. 

Figure ITT1 plots the spatial variation of the ratio of the 
model prediction (case 2) to the observation for the [O i] 
63 fim/ FIR ratio and Z. The two quantities are plotted 
in Fig. El indicating a good correlation. FIR and [O i] 
146 fim are not affected by overlapping because both emis- 
sions are optically thin (r <C 1), and we assumed that the 
results for case 2 can directly be used. The PDR model 
predicts the observed [Oi] 63 fim/ FIR ratio fairly well 
(model/observation ~ 1) for the positions with Z < 1. 
Figure m shows that the regions of d > 3.5 pc have 
Z larger than unity, where the discrepancy between the 
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Fig. 9. Results of the PDR model with [Oi] 63 ^luy/FIR as the input parameter (case 1). a) gas density, b) surface 
temperature, c) model prediction (solid errorbars) and observation (dashed errorbars) of [Oi] 146 fj,m/FIR, and d) 
model prediction (solid errorbars) and observation (dashed errorbars) of [Cll] 158 /an/ FIR. The errors of the model 
prediction include the uncertainties in the observed [Oi] 63 ^m intensities and Go- 



E 

CO 
CD 



— 2 



(a) 



P. 3 



2 



1 

10 



I I 



(c) 



I 1 M hHfl 



800 



cu 600 



400 



E 
=i. 

CO 

in 



O 



200 
14 
12 

1 10 

- 8 

6 
4 
2 




(b) 



(d) 



i s t 



Distance from Be 59 [pc] 
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Fig. 11. Ratio of the model prediction (case 2) to the 
observations for the [O i] 63 /im/FIR ratio (solid errorbars) 
and the overlapping factor Z (dotted errorbars; Eq. (JBJ)). 
The overlapping factors are shifted slightly to the right for 
clarity. 




Fig. 12. Ratio of the model prediction (case 2) to the ob- 
servation of [Oi] 63 /im/FIR against the overlapping fac- 
tor Z. The dashed line indicates the relation calculated by 
Eq. @ with r> 1, while the dotted line indicates that 
shifted by -50% in Z. 



model and the observations is significant (see also Figs.llOb 
and d). The value of Z larger than unity indicates more 
than one clouds in the beam, suggesting that absorption 
of the overlapping clouds is significant for the [O i] 63 /im 
emission. 

Assuming that N PDR clouds, each of which has the 
same optical depth r, overlap on the line of sight, the 
total observed intensity normalized by the intensity from 
a cloud can be given by 

JV-l i _ p -Nr 

3=0 

We assume t > 1 for [Oi] 63 /Ltm. The line of Eq. © 
with t ^> 1 is also plotted in Fig. ^| by substituting N 
by the overlapping factor Z . Figure ITU indicates that the 
correlation can be accounted for by Eq. (J5J if Z is shifted 
by —50% (dotted line). Such a shift can be attributed to 
the uncertainty in the absolute scale of FIR and/or to the 
simple assumption of the model. 

For [Cn] 158 /im, we assume r = 1 and apply Eq. © 
with the shift of —50% in Z to the model results predicted 
from Go and [Oi] 146 /im/ FIR. The discrepancy that the 
model intensity is stronger than the observations for [C n] 
158 /im for d > 3.5 pc, is resolved by taking account of 
the absorption due to the overlapping clouds. The resul- 
tant model prediction for the [Cn] 158 /im emission be- 
comes less than those observed at all the positions. Then 
the excess intensity in the observed [C II] 158 /Ltm over the 
overlapping model prediction can be ascribed to that from 
the ionized gas. FigurelT3lplots the thus estimated fraction 
of the [Cn] 158 /im emission from the ionized gas, indicat- 
ing that a non-negligible fraction (20-70%) of the [C n] 
158 /im emission comes from the ionized gas. The amount 
of the ionized gas contribution to [Cn] 158 /im is within 
a reasonable range expected from the [Nil] 122 /im inten- 



1 

to 

D5 T 
"S 0.8 - T T 

— T 

"d T 

o 

E T T T 

2 0.6 t It 

E I 1 

i . 1 

co 1 T 

in 1 

^ 0.4 - 11 

o III 

o T 

J 0.2 - 1 1 " 

1 , J , 

1 2 3 4 5 

Distance from Be 59 [pc] 

Fig. 13. Estimated fraction of [C n] 158 /im from the ion- 
ized gas. The errorbars include the uncertainties in the 
observed intensities of [Cn] 158 /im and [Oi] 146 /im, and 
Go- 

sity (see Sect. I3.3|) . Thus the simple model of overlapping 
PDR clouds accounts for the observations consistently. 

3.3. Lowly-ionized gas 

[Nil] 122 /im, and part of [Cn] 158 /im and [Sin] 35 /im 
originate from the lowly-ionized gas. Using the observed 
[Nil] 122 /im intensity and the [Cn] 158 /im intensity es- 
timated from the ionized gas in Sect. 13.21 we derive the 
electron density by assuming the abundance ratio of C + 
to N + . Model calculations suggest that in the warm low- 
density photoionized interstellar medium, 81% of nitrogen 
is in N+ and 95% of carbon is in C+ (Sembach et al-HHP. 
The difference of the ionization fractions is smaller than 
that in the estimate of the [C n] 158 /im intensity from the 
ionized gas and we simply assume that the ratio of C + to 
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Fig. 14. Observed [Nil] 122 fim intensity versus estimated 
[Cn] 158 fim intensity from the ionized gas. The dashed 
lines indicate the calculations for the solar abundance gas 



70, 30, and 10 cm 3 (from the upper to lower 



lines) . The dotted lines show those for the ISM abundance 
gas of n e — 10 cm~ 3 (upper line) and at the low density 
limit (lower line) . The electron temperature is assumed to 
be 10 4 K for all cases. 
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Fig. 15. Observed [Nil] 122 intensity versus [Sill] 
35 /im intensity. The upper line indicates the gas of the 
Si abundance of 40% of solar and the low density limit. 
The middle line indicates the gas of the Si abundance of 
30% of solar and n e = 30 cm~ 3 . The lower line indicates 
the gas of the Si abundance of 20% of solar and n e = 70 
cm~ 3 . The electron temperature is assumed to be 10 4 K 
for all the cases. 



N + is equal to the gas phase abundance ratio of C to N. 
FigurelHIplots the observed [N n] 122 fim intensity against 
the [C II] 158 fim intensity from the ionized gas. We assume 
two cases for the abundance, the latest solar abundance 
(log(N/H) = -4.07 and log(C/H) = -3.41; Holweger 
l2TM)l and the ISM abundance (log(N/H) = -4.10 and 
log(C/H) = -3.86; Savage & Sembach IT^jl . Figure d 
suggests that the observed [N n] 122 /im and [C n] 158 /im 
intensities can arise from the ionized gas with the elec- 
tron density of about 30^20 cm -3 for the solar abundance 
and that of less than 10 cm -3 for the ISM abundance for 
T e = 10 4 K. COBE observations suggest that the electron 
density of the [Nil] 122 /im emitting gas in the Galactic 
plane is typically 10-30 cm -3 (Wright et al. I1991fl . which 
is in agreement with the present estimate. We estimate 
the path length of the N + gas in the same manner as for 
the ++ gas. With the strongest [N n] 122 /im intensity 
of 2.9 x lO^Wm^sr 1 and assuming T e = 10 4 K and 
the solar abundance, we derive the length L = 4.2 and 
0.9 pc for n e = 30 and 70 cm -3 , respectively. Because of 
the large uncertainty in n e due to the uncertainty in the 
abundance, it is difficult to conclude whether the lowly- 
ionized gas is distributed in a thin sheet or extending over 
the S171 region. 

The observed [Nil] 122 /xm and [Sin] 35 fim intensities 
give us an estimate of the gas phase silicon abundance. 
The N abundance is assumed to be solar without deple- 
tion and the singly ionized ion fraction is assumed to be 
the same for N and Si. Both lines are also assumed to orig- 
inate from the same ionized region. Silicon and nitrogen 
have quite different ionization potential energies. However, 
model calculations indicate that 81% of nitrogen is in N + 
and 91% of silicon is in Si + in the warm ionized medium 



(WIM) owing to the charge exchange effect (Sembach et 
al. 12000(1 . Taking account of the calculated ionization frac- 
tions the required silicon abundance is lowered by about 
10% compared to the assumption of the same ionization 
fraction. Figure IT5l shows the observed intensities and the 
results of the calculation. S171 shows strong silicon emis- 
sion, which requires log (Si/H) = —4.99, 30 ± 10% of solar 
abundance silicon in the gas phase for n e < 70 cm -3 . 

Figure El shows that the [Si n] 35 /im emission is 
well correlated with the [N Ii] 122 /im emission except for 
the two positions where the [Si n] 35 /im emission is the 
strongest. This correlation suggests that the [Sin] 35 fim 
emission originates mainly from the [Nil] 122 /im emit- 
ting gas. Observations of the Carina region by Mizutani 
et al. ( 2000 ) also suggest a good correlation between [Si n] 
35 fim and [Nil] 122 fim, indicating that the [Sin] 35 /im 
emission comes mostly from the ionized gas. The strong 
[Sin] 35 fim emission of the two points, d ~ 4 pc, may be 
attributed to the contribution from the PDR or the high 
electron density gas. 

The [Sin] 35 fim line has a high critical density (2.6 x 
10 5 cm -3 ) for the H° collision and a high exciting tem- 
perature (413 K). Hence, in the neutral region the [Sin] 
35 fim emission indicates the presence of high-temperature 
and high-density gas. In Sect. 13.21 the density of the PDR 
in S171 is estimated to be 100-1500 cm -3 and the temper- 
ature is 200 K-400 K. These properties indicate that the 
[Si n] 35 fim emission from the PDR should not be signifi- 
cant as Fig. E| indicates. Based on the comparison of the 
observed [Si n] 35 fim emission with the PDR model by 
Hollenbach et al. ifjMTj) . we found that we need by about 
12 times larger silicon abundance in the gas than that as- 
sumed in the PDR model (about 2% of solar abundance) if 
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the strongest [Si n] 35 pm emission (d ~ 4 pc) comes only 
from the PDR. The [Si n] 35 pm emission should also be 
attenuated by overlapping clouds because of its large opti- 
cal depth. Taking account of this effect, we estimate that 
silicon of 65 ± 20% of solar abundance must be in the gas 
phase in the PDR at ~ 4 pc. If a half of the [Si n] 35 /im 
emission at ~ 4 pc comes from the ionized gas and the 
other half from the PDR, then both regions need gaseous 
silicon of about 30% of solar abundance. Therefore the 
large Si abundance in the gas in the S171 region is a se- 
cure conclusion. 

Intense [Si n] emission has been reported in several H n 
regions (Haas et al. [1991), in the Galactic center region 
(Stolovy et al. I1995f) . and in the Carina region (Mizutani 
et al. I2000J) . Silicon and carbon are major constituents 
of interstellar grains (Mathis I199()|l . The gas phase Si 
abundance is about 5% of solar in cool clouds (Savage & 
Sembach ll996fl and the depleted atoms are thought to re- 
side in interstellar grains. The suggested large abundance 
of gaseous silicon in PDRs or Hll regions indicates ef- 
ficient grain destruction processes taking place in these 
regions. Observations of the massive star-forming region 
the Carina nebula indicate the ratio of [Si n] 35 pm to 
[Nil] 122 /jm of about 8 (Mizutani et al. 2000), while it is 
about 4 in S171. In both regions, the estimated electron 
density is in a similar range (10-50 cm -3 ) and thus the 
difference in the ratio can be attributed to the difference 
in the relative abundance. Observations of the starburst 
galaxy NGC 253 also indicate that about a half of the solar 
abundance Si should be in the gas (Carral et al. 1994). The 
present result indicates a slightly smaller value of ~ 30%. 
The degree of destruction may indicate the difference in 
the activity of the regions. Theoretical investigations sug- 
gest that dust grains are destroyed quite efficiently by su- 
pernova shocks in a time scale of 2-4 x 10 s yr (Jones 
et al. IT39H RMfy . Observations by the Goddard High 
Resolution Spectrograph (GHRS) on board the Hubble 
Space Telescope have indicated that Si returns to the gas 
phase relatively quickly, while Fe grains seem to be more 
resistive (Fitzpatrick 199(2 H997f) . The present observa- 
tions support these results and indicate that Si in dust 
grains is easily returned to the gas phase. This is in con- 
trast to the constant gas abundance suggested for oxy- 
gen and carbon (Cardelli et al. 119961 Sofia et al. 119971 
Cartledge et al. 2001 ). The recycling as well as the compo- 
sition of dust grains in interstellar space should be reexam- 
ined based on these recent results of the gas abundance in 
various phases (cf. Tielens lHMl Jones HUM OnakaEH3) ■ 

4. Discussion 

4.1. The absorption of [0\] 63fim and [Cu] 158 pirn 

In the S171 region, strong [Oi] 146 pm emission is ob- 
served compared to [O i] 63 /im at ~ 4 pc, where the ratio 
of [O i] 63/im/146/im is about 5. Small ratios compared 
to the PDR model prediction have been reported also in 
other objects. Liseau et al. ©99) showed that the ratio 



becomes 1-5 in the p Oph cloud and were not able to iden- 
tify any convincing mechanisms to explain the discrep- 
ancy. They took account of the collisional excitation at low 
temperatures and examined the collision coefficients of O 
in detail. The model prediction based on the newly-derived 
collision parameters, however, did not make a significant 
difference from the older one. Thus the uncertainty in the 
collision coefficients is not a major factor for the low [O i] 
63 ^m to 146 /im ratio. The 3 Pi and 3 Po levels in Oi can 
be populated by cascades from higher states that have 
been populated from 3 P2 through the absorption of the 
UV interstellar radiation (Keenan et al. I1994f) . With the 
oscillator strengths given by Morton l|1991fl we calculated 
the detailed balance of O I level populations including this 
effect. The UV pumping is effective at low hydrogen den- 
sities (e.g. n H < 1000 cm" 3 and T < 1000 K) and the [O i] 
63 /im to 146 pm ratio becomes smaller than 15 for temper- 
atures of several hundred K. Taking account of the possi- 
ble optical depth effect, this process can account partly for 
the obtained low ratio. However, the low ratio is observed 
at the surface of the molecular region, while the UV pump- 
ing is expected to be effective in the vicinity of the star. 
The spatial distribution of the line ratio suggests that the 
UV pumping is not a major mechanism for the low ratio. 
As described in Sect. 13.11 dust extinction does not affect 
the observed line intensity for A > 50 /im and thus does 
not affect the observed [Oi] line ratio. Caux et al. ^ 1999 ) 
claimed that the low value of [O i] 63 /im to 146 /im in the 
p Oph cloud is due to the presence of a very large column 
density of atomic oxygen, which makes the 63 /im line op- 
tically thick. Giannini et al. (2000) showed that the com- 
plex site of massive star formation, NGC 2024 (Orion B, 
W12), has the [Oi] 63 pm to 146 pm ratio of about 5 and 
suggested that the [O i] 63 pm line is strongly absorbed 
by the cold foreground gas. Absorption of the [Oi] 63 pm 
line has been reported directly by high spectral resolution 
KAO observations toward star-forming regions (Poglitsch 
et al. lTMol Kraemer et al. HM5|l . and by ISO/LWS obser- 
vations toward Sagittarius B2 (Baluteau et al. 1997) and 
W49N (Vastel et al. l27TO . 

Carbon atoms are in CO in cold molecular clouds, and 
thus the [C n] 158 pm emission will not be absorbed in cool 
clouds as efficiently as [Oi] 63 pm. However, [Cu] 158 pm 
is marginally optically thick, r ~ 1, in the PDR (Kaufman 
et al. lT§59T> . Boreiko & Betz (P§5) observed NGC 6334 
at high spectral resolution with the KAO and indicated a 
self- absorption profile of the [C n] 158 pm emission at sev- 
eral positions. NGC 6334 is a star- forming region where 
self- absorption of [O i] 63 pm is also observed (Kraemer 
et al. UnnSJ). Boreiko & Betz ifTM?}) further directly de- 
tected the absorption feature of [C II] 158 pm in the star- 
forming region, NGC 3576, and suggested that the [Cu] 
158 pm emission is absorbed by a foreground cloud. The 
present analysis indicates that the model of overlapping 
PDR clouds accounts for observations sufficiently well. 
The resultant [C il] 158 pm emission from the ionized gas 
is in the reasonable range, supporting the present simple 
model. A detailed study of the p Oph region, where the 
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contribution to [C n] 158 /im from the ionized gas is neg- 
ligible and thus the observed line intensities can be more 
directly compared to model predictions, also suggests that 
the [Ol] 63 (j.m, 146 ^m, and [Cn] 158 /im intensities are 
able to be accounted for consistently by the present model 
with the overlapping factor (Okada et al. 2003). 

These results suggest that the absorption of [O i] 63 /im 
and [C n] 158 /im in clouds on the line of sight is the most 
likely cause for the difference between the PDR model 
predictions and the present observations. The overlap- 
ping model does not necessarily assume separate clumpy 
clouds, but can also be applied if a thick cloud is viewed in 
an edge-on configuration. In the PDR model by Kaufman 
et al. l|1999f) . the turbulent velocity dispersion is assumed 
to be 1.5kms _1 , which corresponds to the line width of 
2.5 km s -1 in FWHM. In S171, the 13 CO emission seems 
to have two components, one with the width of 1.0 km s" 1 
and the other with SAkms- 1 width (Yang & Fukui UM^I . 
Thus if the clouds are overlapping, self-absorption in [O i] 
63 /im and [Cn] 158 /im can be expected to some extent 
in SI 71. The validity of the present model of overlap- 
ping PDRs can be examined in detail by high spectral- 
resolution observations. The present analysis suggests that 
the best diagnostic line for PDRs is [O i] 146 /im, and that 
FIR is an important parameter to estimate the degree of 
the overlapping on the line of sight. 

4.2. Comparison with observations of other H II 
regions 

Observations of compact Hn regions suggest that 
their electron density is in the range 200-10000 cm -3 
(Moorwood et al. ITMHI Watson et al. HMT1 Martin- 
Hernandez et al-IHlSl- Observations of G29.96-0.02 sug- 
gest the presence of a dense (n c ~ 57000 cm~ 3 ) core and a 
diffuse component (n c ~ 680 cm' 3 ; Morisset et al. '2002). 
The H II region- molecular cloud complex of NGC 2024 has 
been shown to have n c — 1200-1500 cm -3 (Giannini et 
al. |2000). The present results of S171 indicate the lower 
end of the electron density in the range of these observa- 
tions, suggesting the presence of the diffuse H n region in 
S171. This is supported by the extended distribution of 
the [O in] and [N in] lines over several pc. Observations of 
the Carina nebula also indicate the presence of the low- 
density (< 100 cm~ 3 ) highly-ionized gas extending over 
several tens pc (Mizutani et al. I2THJ2T). The presence of 
low-density highly-ionized gas in star-forming regions may 
not be an uncommon phenomenon. 

In the Carina nebula, however, the [O i] 63 /im emis- 
sion is weak and is about one third of the [C il] 158 /im 
emission. The [O i] 146 /im emission is also weak and has 
been detected only at several positions (Mizutani et al. 
2000). Compared to the Carina nebula, SI 71 has a large 
neutral gas density. On the other hand, observations of 
NGC 2024 suggest a much higher density (~5x 10 5 -10 6 
cm" 3 ) for the PDR (Giannini et al-EEE!)- S171 may be at 



the middle stage in between dense and young star-forming 
regions and evolved diffuse PDRs. 

The fraction of the [C n] emission from the ionized gas 
is suggested to be less than 30% both in the NGC 2024 
and the Carina regions. The present observations suggest 
a higher ratio of 20-70% of the [C n] emission that comes 
from the ionized gas. This ratio has a large uncertainty 
because it depends on the overlapping model, but if real, 
the large fraction may be partly due to the geometrical 
effect: in SI 71 the ionized gas and PDR overlap on the line 
of sight with a large degree. Observations of the gas motion 
should provide useful information to further investigate 
the geometry of the S171 region. 

5. Summary 

We observed S171 with the LWS, the SWS, and the PHT- 
S on board ISO. We have detected 7 forbidden lines, 
[Om] 52 /im, [Nm] 57 /im, [Ol] 63 /im, [Om] 88 /im, [Nil] 
122 /tm, [Ol] 146 /im, and [Cn] 158 /im, and far-infrared 
continuum with the LWS, [Sill] 35 /im and H2 9.66 /im 
with the SWS. The PHT-S observations gave only upper 
limits on the H2 9.66 /im and 6.9 /im lines. We attribute 
the observed lines to three phases - highly-ionized gas, 
lowly-ionized gas, and PDR gas - and derived the physi- 
cal properties of each phase. 

For the highly-ionized gas, we have estimated the elec- 
tron density from the [Om] lines as n e < 150 cm" 3 . The 
path length of the ++ gas along the line of sight is es- 
timated to be 9 x (20cm~ 3 /n e ) 2 pc, indicating that the 
highly-ionized gas is associated with S171. The abundance 
ratio of N ++ /0 ++ is 0.2 ± 0.1, in agreement with the ob- 
served ratio of Hn regions at the galactocentric distance 
of about 8 kpc. 

The [Ol] 63 /im, 146 /im and [Cn] 158 /im line intensi- 
ties around the interface region between the ionized and 
molecular gases cannot be accounted for consistently by 
the current PDR model. The [0 1] 146 /im emission is very 
strong, and the [0 1] 63 fim to 146 fim ratio is too small. 
The [Cn] 158 /im emission also seems to be too weak. 
We conclude that the model based on the observed [0 1] 
146 /im/ FIR ratio gives the most consistent results with 
the other observations. We propose a model in which the 
[0 1] 63 /im and [C 11] 158 /im emissions are attenuated by 
overlapping PDR clouds. We introduce the overlapping 
factor Z defined as the ratio of the observed FIR to the 
far-infrared intensity of a PDR cloud estimated from the 
temperature of the FIR continuum emission. The observed 
intensities of the three lines can be accounted for by the 
simple model of overlapping PDRs. We obtain the gas den- 
sity n ~ 100-1500 cm -3 and the temperature T ~ 200- 
400 K. The present analysis indicates that the [0 1] 63 /im 
line emission must be used with caution in diagnostics of 
PDRs and the [Ol] 146 /im line may be the best parame- 
ter for the diagnostics. A comparison of the observations 
with the model prediction indicates that 20-70% of the 
[C 11] 158 /im comes from the ionized gas. 
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We have obtained the electron density n e = 30I 20 
cm~ 3 for the solar abundance or less than 10 cm~ 3 for 
the ISM abundance for the lowly- ionized gas. [Sin] 35/im 
is quite strong and the observed ratio of [Sin] 35/im to 
[Nil] 122 fim indicates the silicon abundance in the diffuse 
ionized gas to be about 30% of solar. This value is much 
larger than the typical value in interstellar diffuse clouds 
(~ 5% of solar), suggesting the efficient dust destruction 
in the ionized region. The present observations indicate 
the importance of the [Si n] line in future observations to 
investigate the dust destruction and the circulation pro- 
cesses of the interstellar medium. 
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